Dynamical Quantum Phase Transitions in the Transverse Field Ising Model 



M. Heyl 

Department of Physics, Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, 
Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, 80333 Munich, Germany and 
Institut fur Theoretische Physik, Technische Universitat Dresden, 01062 Dresden, Germany 

A. Polkovnikov 

Department of Physics, Boston University, 590 Commonwealth Ave., Boston, MA 02215, USA 

S. Kchrein 

Department of Physics, Georg- August- Universitat Gottingen, 
Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 

A phase transition indicates a sudden change in the properties of a large system. For temperature- 
driven phase transitions this is related to non-analytic behavior of the free energy density at the 
critical temperature: The knowledge of the free energy density in one phase is insufficient to predict 
the properties of the other phase. In this paper we show that a close analogue of this behavior can 
occur in the real time evolution of quantum systems, namely non-analytic behavior at a critical time. 
We denote such behavior a dynamical phase transition and explore its properties in the transverse 
field Ising model. Specifically, we show that the equilibrium quantum phase transition and the 
dynamical phase transition in this model are intimately related. 



Phase transitions are one of the most remarkable phe- 
nomena occurring in many-particle systems. At a phase 
transition a system undergoes a non- analytic change of 
its properties, for example the density at a tempera- 
ture driven liquid-gas transition, or the magnetization 
at a paramagnet-ferromagnet transition. What makes 
the theory of such equilibrium phase transitions particu- 
larly fascinating is the observation that a perfectly well- 
behaved microscopic Hamiltonian without any singular 
interactions can lead to non-analytic behavior in the ther- 
modynamic limit of the many-particle system. In fact, 
the occurrence of equilibrium phase transitions was ini- 
tially a puzzling problem because one can easily verify 
no go theorems for finite systems, therefore the thermo- 
dynamic limit is essential pQ. 

Today the theory of equilibrium phase transitions is 
well established, especially for classical systems under- 
going continuous transitions, where the powerful tool of 
renormalization theory bridges the gap from microscopic 
Hamiltonian to universal macroscopic behavior. On the 
other hand, the behavior of non-equilibrium quantum 
many-body systems is by far less well understood. Recent 
experimental advances have triggered a lot of activity in 
this field [3] , like the experiments on the real time evolu- 
tion of essentially closed quantum systems in cold atomic 
gases [31 Q]. The experimental setup is typically a quan- 
tum quench, that is a sudden change of some parameter 
in the Hamiltonian. Therefore the system is initially pre- 
pared in a non-thermal superposition of the eigenstates 
of the Hamiltonian which drives its time evolution. 

From a formal point of view, there is a very suggestive 
similarity between the canonical partition function of an 
equilibrium system 

Z(0)=Tre-^ H (1) 



and the overlap amplitude of some time-evolved initial 
quantum state l^) with itself 

G(t) = (^e-^i) (2) 

This leads to the question whether some analogue of tem- 
perature (/3)-driven equilibrium phase transitions in 
exists in real time evolution problems. In the theory of 
equilibrium phase transitions it is well established that 
the breakdown of the high-temperature (small /3) ex- 
pansion indicates a temperature-driven phase transition. 
Likewise, we propose the term dynamical phase transition 
for non-analytic behavior in time, that is the breakdown 
of a short time expansion in the thermodynamic limit 
at a critical time. In this paper we study this notion of 
dynamical phase transition in the one dimensional trans- 
verse field Ising model, which serves as a paradigm for 
one dimensional quantum phase transitions [5]. It can 
be solved exactly, which permits us to establish the exis- 
tence of dynamical phase transitions that are intimately 
related to the equilibrium quantum phase transition in 
this model. 

Our key quantity of interest is the boundary partition 
function 

Z{z) = {ne- zH m (3) 

in the complex plane z € C. For imaginary z = it this 
just describes the overlap amplitude pi). For real z = R 
it can be interpreted as the partition function of the 
field theory described by H with boundaries described by 
boundary states separated by R [BJ. In the thermo- 
dynamic limit one defines the free energy density (apart 
from a different normalization) 

/(*) - - Hm 1 In Z{z) (4) 
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where N is the number of degrees of freedom. Now sub- 
ject to a few technical conditions [1] one can show that 
for finite N the partition function ^ is an entire func- 
tion of z since inserting an cigenbasis of H yields sums 
of terms e~ zEj , which are entire functions of z. Accord- 
ing to the Weierstrass factorization theorem [7] an entire 



function with zeroes 



Za G 



can be written as 



^)=e^n(i-£) 



(5) 



with an entire function h(z). Thus 



lim — 



fc(*) + X> (1 



(6) 



and the non-analytic part of the free energy density is 
solely determined by the zeroes Zj . A similar observation 
was originally made by M. E. Fisher [T], who pointed 
out that the partition function is an entire function 
in the complex temperature plane. This observation is 
analogous to the Lee- Yang analysis of equilibrium phase 
transitions in the complex magnetic field plane [8]. For 
example in the 2d Ising model the Fisher zeroes in the 
complex temperature plane approach the real axis at the 
critical temperature z = /3 C in the thermodynamic limit, 
indicating its phase transition [5]. 

We now work out these analytic properties explicitly 
for the one dimensional transverse field Ising model (with 
periodic boundary conditions) 



H(g) 



JV-l 

E* 2 

»=i 



'i+l 



(7) 
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For magnetic field g < 1 the system is ferromagnetically 
ordered at zero temperature, and a paramagnet for g > 1 
[5] . These two phases are separated by a quantum critical 
point at g = g c = 1. The Hamiltonian Q can be mapped 
to a quadratic fermionic model [TUHT2"] 



N-l 



N 

-gJ2 c U (8) 

»=i 



Diagonalization yields the dispersion relation e k (g) = 

\J (g — cos k) 2 + sin 2 k. 

In a quantum quench experiment the system is pre- 
pared in the ground state for parameter g , — 
\^Gs(go)}: while its time evolution is driven with a 
Hamiltonian H{g\) with a different parameter g\. In the 
sequel we will first analyze quench experiments in the set- 
ting of the fermionic model Q. A subtle difference occurs 
when thinking in terms of the spin model Q since in the 
ferromagnetic phase the ground state of the spin model 
is twofold degenerate, while the fermionic model always 
has a unique ground state. We will say more about this 
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Figure 1: Lines of Fisher zeroes for a quench within the same 
phase go = 0.4 — ^ g\ = 0.8 (left) and across the quantum 
critical point go = 0.4 — > g 1 — 1.3 (right). Notice that 
the Fisher zeroes cut the time axis for the quench across the 
quantum critical point, giving rise to non-analytic behavior 
at t*n (the times t* are marked with dots in the plot). 



later. Taking the ground state of the fermionic model in 
Eq. ([8| as the initial state |^) the free energy density Q 
describing this sudden quench g — > g 1 can be calculated 
analytically [J] yielding 



fgo,gi ( z ) 



* dk 
2^ 



In (cos 2 (t> k + sin 2 cj> k e ~ 2zek ^ 



(9) 

Here 4> k = k (g ) - k (gi), and tan(26 k (g)) = smk/[g - 
cos/c) , 0k(g) G [0,7r/2] . In ((9j) we have ignored an un- 
interesting additive contribution z EGs(gi)/N that de- 
pends on the ground state energy of H{gi). 



In the thermodynamic limit the zeroes of the partition 
function in the complex plane coalesce to a family of lines 
labeled by a number n G Z 



z n (k) 



1 



2e/c (51) 



(lntan 2 fe +in(2n+ 1)) (10) 



»fc=0 



The limiting infrared and ultraviolet behavior of the 
Boboliubov angles 

quench in same phase 
7r/4 quench to/from quantum critical point 
7r/2 quench across quantum critical point 

aW = (11) 

immediately shows that the lines of Fisher zeroes cut 
the time axis for a quench across the quantum crit- 
ical point (Fig. [T]) since then limfe_j.oRe z n (k) = 00, 
lim^^Tr Re z n (k) = —00. In fact, the limiting behav- 



ior (11) remains unchanged for general ramping proto- 



cols |14) . 
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The free energy density Q is just the rate function 
of the return amplitude G(t) = exp[— N f{it)]. Like- 
wise for the return probability (Loschmidt echo) L(t) = f 
\G{t)\ 2 = exp(-iV7(t)) one has l(t) = f(it) + f(-it). 
The behavior of the Fisher zeroes for quenches across 
the quantum critical point therefore translates into non- 
analytic behavior of the rate functions for return am- 
plitude and probability at certain times t* . For sudden 
quenches one can work out these times easily 

t* n =t* (n+^j , n = 0,1,2,... (12) 

with f* = 7r/efc«(<7i) and k* determined by cosfc* = (1 + 
9o 9i)/(9o+9i)- We conclude that for any quench across 
the quantum critical point the short time expansion for 
the rate function of the return amplitude and probability 
breaks down in the thermodynamic limit, analogous to 
the breakdown of the high-temperature expansion at an 
equilibrium phase transition. In fact, the non-analytic 
behavior of l(t) at the times t n has already been derived 
by Pollmann et al. [15j for slow ramping across the quan- 
tum critical point. For a slow ramping protocol £fc*(gi) 
becomes the mass gap m(gi) = \gi — 1| of the final Hamil- 
tonian, but in general it is a new energy scale generated 
by the quench and depending on the ramping protocol. 
In the universal limit for a quench across but very close 
to the quantum critical point, g± = I + 5, \S\ <C 1 and 
fixed go, one finds e^* {gi)/m(gi) <x 1/ y/\5\. Hence in this 
limit the non-equilibrium energy scale e^* becomes very 
different from the mass gap, which is the only equilibrium 
energy scale of the final Hamiltonian. 

The interpretation of the mode k* follows from the ob- 
servation n(k*) — 1/2, where n{k) is the occupation of 
the excited state in the momentum fc-mode in the eigen- 
basis of the final Hamiltonian Hf(gi). Modes k > k* 
have thermal occupation n(k) < 1/2, while modes k < k* 
have inverted population n{k) > 1/2 and therefore for- 
mally negative effective temperature. The mode k* cor- 
responds to infinite temperature. In fact, the existence 
of this infinite temperature mode and thus of the Fisher 
zeroes cutting the time axis periodically is guaranteed for 
arbitrary ramping protocols across the quantum critical 
point. For example, for slow ramping across the quantum 
critical point the existence of this mode and the negative 
temperature region in relation to spatial correlations was 
discussed in Ref. [To] . 



One measurable quantity in which the non-analytic be- 
havior generated by the Fisher zeroes appears naturally 
is the work distribution function of a double quench ex- 
periment: We prepare the system in the ground state of 
H(go), then quench to H(gi) at time t = 0, and then 
quench back to H(go) at time t. The amount of work W 
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Figure 2: The bottom plot shows the work distribution func- 
tion r(w,t) for a double quench across the quantum critical 
point (go = 0.5, gi — 2.0). The dashed line depicts the ex- 
pectation value of the performed work, r(w, t) = 0. The top 
plot shows various cuts for fixed values of the work density w. 
The line to = is just the Loschmidt echo: Its non-analytic 
behavior at t* n becomes smooth for w > 0, but traces of the 
non-analytic behavior extend into the work density plane. In 
this respect work density plays a similar role to temperature 
in the phase diagram of an equilibrium quantum phase tran- 
sition. 

performed follows from the distribution function 

W*)=E^-(^f-W)) K^l*i(*))l a 

(13) 

where the sum runs over all eigenstates \Ej) of the ini- 
tial Hamiltonian H(g ). It obeys a large deviation form 
P(W,t) ~ e -Nr(.w,t) with a rate f unct i 0n r ( Wt t) > de- 
pending on the work density w — W/N. In the thermo- 
dynamic limit one can derive an exact result for r(w,t): 
According to the Gartner-Ellis theorem [17] it is just the 
Legendre transform 

- r(w, t) = inf (wR - c(R, t)) (14) 

_R.GR 

where 

fx (Ik f 

c{R,t) = - J — ln(l + sin 2 (2(/) fc )sm 2 (e fe ( 5l )t) 

x ( e - 2e »<»0 B -l)) (15) 

is the rate function for the cumulant generating func- 
tion of the work distribution function, C(R,t) = 
J dW P{W 7 t)e- RW = e - N <*t). In Fig. Fig. |] we 
show r(w,t) for a quench across the quantum critical 
point. For w — it just gives the return probability to 
the ground state, r(w = 0, t) = lit), therefore the non- 
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analytic behavior at the Fisher zeroes shows up as non- 
analytic behavior in the work distribution function. How- 
ever, from Fig. [2] one can see that these non-analyticities 
at w = also dominate the behavior for w > at t* , cor- 
responding to more likely values of the performed work. 
The suggestive similarity to the phase diagram of a quan- 
tum critical point, with temperature being replaced by 
the work density w, motivates us to call this behavior 
dynamical quantum phase transitions. Notice that ex- 
perimentally the work density can be lowered by post- 
selection [T5] . 

So far we have analyzed the quench dynamics in terms 
of the fermionic model ph. When thinking in terms of the 
transverse field Ising model Q, all results carry over for 
quenches starting in the paramagnetic phase since then 
the spin ground state is unique. Specifically, one finds 
the non-analytic behavior in the Loschmidt echo and the 
work distribution function for quenches from the param- 
agnetic to the ferromagnetic phase. For quenches orig- 
inating in the ferromagnetic phase, the Loschmidt echo 
calculated above corresponds to working in the Neveu- 
Schwarz sector [21j . which amounts to an unphysical su- 
perposition of spin up and spin down ground states in 
the spin language. However, looking at the experimen- 
tally relevant quantity work distribution function, one 
derives the same result in the thermodynamic limit as 
above when starting from either of the two degenerate 
ferromagnetic ground states. Specifically, one obtains the 
non-analytic behavior in P(w — 0, t) at the critical times 
(12 1 for quenches from the ferromagnetic to the param- 




time (t — tj)/t* 

Figure 3: Dynamics of the magnetization after the quench. 
The bottom plot shows the longitudinal magnetization for 
various quenches across the quantum critical point. The time 
axis is shifted by a fit parameter t v and one can see that the 
period of the oscillations is the time scale t* ( 12 1. The upper 
plots show the magnetization dynamics in the y — z-plane for a 
quench across the quantum critical point go — 0.3 —> g\ = 1.4 
(left) and a quench in the ordered phase go = 0.3 — > gi = 0.8 
(right). For better visibility the magnetization is normalized 

,,{t)/s/sl{t)+ s%{t). Notice the 



to unit length: s y 



, (t) = Sy 



agnetic phase 



Larmor precession for the quench across the quantum critical 
point, while the dynamics for the quench in the ordered phase 
is asymptotically just an exponential decay [19] . 



Interestingly, the non-equilibrium time scale ( 12 ) also 



plays a role in the dynamics of a local observable after 
the quench. We have calculated the longitudinal mag- 
netization by numerical evaluation of Pfafnans [1]. For 
quenches within the ordered phase it is known analyti- 
cally [ini HO] that the order parameter decays exponen- 
tially as a function of time, which is expected since in 
equilibrium one only finds long range order at zero tem- 
perature (g < 1). For a quench across the quantum crit- 
ical point an additional oscillatory behavior is superim- 
posed on this exponential decay, see Fig. [3j Notice that 
the behavior of the magnetization remains perfectly an- 
alytic, but the period of its oscillations agrees exactly 
(within numerical accuracy) with the period t* of Fisher 
times. A conjecture consistent with our observation was 
also formulated in Ref. 21j. A better understanding of 
this observation will be the topic of future work. At low 
energies the oscillatory decay transforms into real-time 
nonanalyticities at the Fisher times using the concept of 
post-selection allowing to observe the dynamical phase 
transitions in local observables |13j . 

Summing up, we have shown that ramping across 
the quantum critical point of the transverse field Ising 



model generates periodic non-analytic behavior at cer- 
tain times t* n . This breakdown of the short time expan- 
sion is reminiscent of the breakdown of a high tempera- 
ture expansion for the free energy at an equilibrium phase 
transition. We have therefore denoted this behavior dy- 
namical phase transition. Notice that there are other re- 
lated but not identical notions of dynamical phase tran- 
sitions, for example a sudden change of the dynamical 
behavior of an observable as a function of some control 
parameter [22, 23 , or qualitative changes in the ensem- 
ble of trajectories as a function of the conjugate field of 
a dynamical order parameter |24j . 

For quenches within the same phase (including to/from 
the quantum critical point) the lines of Fisher zeroes lie 
in the negative half plane, Re Zj(k) < (Fig. IT]) . Hence 
the knowledge of the equilibrium free energy f(R) on the 
positive real axis completely determines the time evolu- 
tion by a simple Wick rotation. This is no longer true for 
a quench/ramping protocol across the quantum critical 
point since then the lines of Fisher zeroes cut the complex 
plane into disconnected stripes, Fig.[TJ Knowing f(R) for 
R > does not determine the time evolution for t > t$. 
In this sense non-equilibrium time evolution is no longer 
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described by equilibrium properties. 
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Loschmidt matrix 

For Hamiltonians with symmetry-broken ground states 
the definition of the return amplitude in Eq. ^ in the 
main text shares an ambiguity: there is not a unique 
ground state and therefore initial state for the 

nonequilibrium time evolution. Moreover, the return am- 
plitude turns out to depend crucially on the precise choice 
of initial state in the degenerate ground state manifold 
for quenches across the quantum critical point in the 
Ising model for go < 1 as will be shown in detail be- 
low. This is of particular importance because the ground 
state \^Gs(go)) of the fermionic version of the Ising chain 
in Eq. ^ is not identical to one of the symmetry bro- 
ken ferromagnetic ground states but rather to some su- 
perposition. The aim of this supplementary material is 
to demonstrate that, although the return amplitude de- 
pends on the precise choice of initial state, the experi- 
mentally relevant quantity, i.e., the work distribution of 
the double quench in Eq. (13), does not. As a conse- 



quence, the free energy density of the fermionic model 
in Eq. ^ is the relevant quantity although \^gs{9o)) is 
not the initial state of the quantum quench protocol. 

In case of symmetry broken ground states we propose 
the following generalization for the return amplitude that 
we will term the Loschmidt matrix: 



C{t) = 



($ 1 | e -"'/ t |tf 1 ) (*i|e-*' H / t |*2) 



(16) 



shown here for the case of two symmetry broken ground 
states l^i) and | \P 2) - The generalization to more than 
two states is straightforward. This matrix depends on the 
particular choice of the basis wave functions \^\) and 
1^2) m the degenerate ground state manifold. From a 
physical point of view, there is, however, a unique choice: 
the symmetry broken ground states. In case of vanishing 
transverse magnetic field in the Ising model this will be 
the fully polarized states = | t • • • t) or |\& 2 ) = | -I 
... I), for example. 

Note that neither nor ^2) are the ground 

state \^Gs(go)) °f the fermionized model in Eq. ^ 
in the ferromagnetic phase. This can be seen di- 
rectly by noticing that l^csido)} carries no magnetiza- 
tion (y GS (g )\af |*gs(.9o)> = 0. Concluding \V GS (g )) 
rather is a superposition of and \^2)- Fortu- 

nately, for most of the quantities considered in the 
literature such as the spin-spin correlation function 
crfa^ this makes no difference. This actually al- 
lows for the calculation of the absolute value of the 
magnetization even in the fermionic language using 
the cluster decomposition (4 , i/2(ffo)|cf |^ / i/2(5 r o)) 2 = 



G 



lim m _ >00 (* G5 (5o)|crf +m crf |*gs(#o)) J]. Some quanti- 
ties, however, such as the magnetization itself and the 
return amplitudes (as will be shown in detail below) the 
precise choice of the initial state can be crucial. 

In the following we demonstrate the dependence of the 
return amplitude on the initial state for quenches from 
the ferromagnetic to the paramagnetic phase. We illus- 
trate this subtle behavior for the particular example of a 
quench from go = to g\ >■ 1. This case can be solved 
exactly and directly for the initial spin Hamiltonian in 
Eq. ([7]) without using the mapping onto the fermionic 
language. For such a quench the Loschmidt matrix is 
given by 



with 



C(t) = 



/n(i) 
/12W 



-Nfn(t) e -Nf 12 (t) 
-Nhi{t) P -Nf 22 (t) 



/22(t) = ]og(|cos(sif/2)|) 
/ 2 i(i)=log(|sin( ffl t/2)|) 



(17) 



(18) 



Here, we have ignored imaginary parts as they are not im- 
portant for the discussion below. Note the large deviation 
scaling of the Loschmidt matrix elements that are expo- 
nentially suppressed with system size N. If Re[fu(t)] > 
Re[/i2(t)] the Loschmidt matrix is effectively diagonal in 
the thermodynamic limit N ^> 1 whereas in the oppo- 
site case it only carries off-diagonal entries. For small 
times the diagonal elements typically dominate. But for 
sufficiently large times the off-diagonal components can 
build up indicating the possibility of transitions between 
the two different symmetry broken ground states. As 
we explain below this switching underlies the dynamical 
phase transition we are discussing. 

The return amplitude G{t) for the ground state 
|*Gs(<7o)) of the fermionic model (which is a superpo- 
sition of |*i) and 1^2)) and thus the free energy density 
f(t) in Eq. ^ is given by 



G(t) = e-^W = (^Gs(9o)\e- ltH n^Gs(9o)) 



Loschmidt matrix changes from the diagonal to the off- 
diagonal sector and vice versa. For the free energy den- 
sity f(t) this implies the following result 



/(*) 



fn(t) , if t| n _j <t <t* 2n . . 

f 12 (t) ,i{t* 2n <t<t* 2n+1 > n£Z - (20) 



Concluding, the singularity in f(t) is in fact not because 
the diagonal free energies are singular but rather because 
at the Fisher times t* the free energy switches between 
fu(t) and fi2(t). Physically it implies that at the first 
Fisher time, for example, the return probability becomes 
dominated by the transition to a different magnetization 
sector. 

Although the precise choice of initial state can be cru- 
cial for the return amplitude, in the following we demon- 
strate that the work distribution function P(W,t) of a 
double quench defined in Eq. ( 13 1 is independent of all 



these subleties. This is important as P(W, t) is the ex- 
perimentally measurable and therefore relevant quantity. 
In the zero work limit one obtains 



Jm P(W, t) = |<*i |*o (t)> I 2 



l(* 2 |*o(i))| 2 (21) 



with 1*0(0) = e ~ 1 |*o) the time evolved initial state 
|*o)- Irrespective of the precise choice of |*o) m the 
ground state manifold (could also be one of the symme- 
try broken ground states) P(W — > 0,t) always contains 
contributions from both magnetization sectors: 



P(W ->(),*) = e -2JVRc[/n(t)] +e -2JVRo[/ 12 (t)] = 



= e -2JVRe[/(t)]_ 



(22) 



For the derivation of this equality the large deviation 
scaling of the probabilities is essentially important. Al- 
though |*Gs(ffo)) is n °t the correct initial state for the 
quantum quench (it is not one of the symmetry broken 
ground states) its associated return amplitude with rate 
function f(t) fully determines the experimentally rele- 
vant work distribution function: 



e -Nfn(t) +g -JV/i2(t)_ 



(19) 



P(W^0) = |(* G s(3o)|e-^*|*Gs(so))| 



(23) 



Due to the large deviation scaling of the Loschmidt ma- 
trix elements the precise coefficients of the superposition 
are subleading and therefore do not appear in this equa- 
tion. This is the case as long as no coefficient is exactly 
zero or scales exponentially with system size. The diag- 
onal elements of the free energy densities fu(t) = /22(*) 
are analytic at the Fisher times t* n = t*(n + 1/2) with 
t* = 7r/<7i. Instead, they show logarithmic singularities 
at times t = 2t*(n + 1/2) ^ t* n for all n,m <E Z. How- 
ever, at the Fisher times Re[/n(t*)] = Re[/i2(t£)], i- e -' 
sin(|£/it* /2|) = cos(|<?it*/2|), and both the diagonal and 
off-diagonal contributions become identical. This signals 
a critical point where the dominant contribution to the 



The analysis in this supplementary material has illus- 
trated the dependence of the return amplitude on the 
choice of initial state for the extreme quench from go = 
to g 1 3> 1. Although the return amplitude depends on 
this choice, the work distribution for a double quench 
does not. 

It is possible to also include 1/gi corrections in the cal- 
culation of the return amplitude using the original spin 
language for large but finite magnetic fields g\. It turns 
out that this improved treatment fully supports the con- 
clusions drawn above. This analysis, however, is beyond 
the scope of this supplementary material and will be pre- 
sented elsewhere. 
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Postselection of observables 

One route towards the experimental observation of the 
dynamical phase transitions in local observables is the 
idea of post-selection. Let us consider the double-quench 
experiment as discussed in the main text but with a 
slightly additional twist. Namely we assume that the sys- 
tem is prepared in the ground state of the initial Hamil- 
tonian (actually the initial state can be arbitrary) then 
the Hamiltonian of the system is quenched to the final 
Hamiltonian Hf for the time t and then quenched back 
to the initial Hamiltonian Hi. After that we allow the 
system to relax to the diagonal ensemble [2] , which using 
a different language also means measuring the energy of 
the system |3j. In the resulting diagonal state we mea- 
sure the expectation value of an arbitrary observable O 
as a function of Hf and t. Note that for observables 
commuting with Hi like energy the projection to the di- 
agonal ensemble is not affecting the result, while for the 
observables which do not commute with Hi such projec- 
tion makes a difference. 

The expectation value of the observable O is then 

{0{t)) = £IW*)WI>I°I"> = $>„(i)(n|0|n>, 

n n 

(24) 

where \n) denotes the energy eigenstates of Hi and 
Pn(t) = \(ip(t)\n)\ 2 are the probabilities of occupation 
of these states. Within a continuum description the sum 
above can be formally written as a continuous integral 
over energies 

(0(t)) = J dEP(E,t)0(E,t), (25) 

where P(E,t) = J2 n p n (t)5(E — E n ) is the work distri- 
bution function for a double quench and 

°( E ^ = E {nP\n)5{E E n ). (26) 

Postselection. The idea of postselection is that one can 
artificially skew the energy distribution P(E, t) by for ex- 
ample disregarding instances with energy above the cer- 
tain threshold. From computational purposes it is more 
convenient, however, to skew the distribution by multi- 
plying it by an exponential factor exp(— f3E), where /? 
plays the role of postselected temperature. We thus can 
formally define 

Pt(E,t) = J^P(E,t)e-? E , (27) 

where Zs(t) is the postselected partition function: 

Z0) = [ dEP{E,t)e-? E = X>a„(t), (28) 
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Figure 4: (Color online) Dynamics of the post-selected mag- 
netization after a quench from gi = 0.5 to gj = 1.5. While the 
blue (solid) curve shows the order parameter dynamics after a 
quench, the green (dashed) curve shows the zero-temperature 
limit f3 — > oo of the post-selected magnetization. 

with Pp n {t) = Pn{t) exp[—/3E]. Then the post-selected 
expectation value of the observable will read: 

<0(% = 5>&>|0|n>. (29) 

n 

yielding in the continuum description: 

(0(% = J dEPp(E,t)Op(E,t). (30) 

We will not analyze in detail this postselection procedure 
here and its similarities and differences with thermody- 
namics (which is equivalent to preselection in our lan- 
guage) since this lies beyond the scope of our work. We 
only note that as the postselected temperature T = 
is lowered we are effectively projecting the observable 
to the ground state manifold. In Fig. [4] we show the 
postselected magnetization in this low temperature limit 
reavealing jumps between the different symmetry bro- 
ken sectors, which are located precisely at the Fisher 
times. Note that this allows one to observe real-time 
nonanalyticities in local observables. We anticipate that 
in this limit the magnetization dynamics will be deco- 
herence free at the expence of effectively excluding large 
amounts of data from the analysis. This point will be a 
subject of our future work. 

Postselection and the complex time return 
amplitude 

Let us make another brief point where postselection 
can be used to obtain nontrivial results. We now go back 
to a single quench of the Hamiltonian from Hi to Hf. 
Then as it was shown in Ref. [Jj up to the phase factor 
the return amplitude G{t) is the Fourier transform of the 



work distribution after the quench: 

G(t) = J dWP(W)exp[iWt], (31) 

where 

P(W) = J2 m\n f )\ 2 S(E r { -E - W), (32) 

n 

where \n^) are the energy eigenstates of the final Hamil- 
tonian and E^ n are the corresponding eigenenergies. 

Using similar considerations as in the previous section 
we can define the postselected work distribution function 
(but this time with respect to the eigenstates of the final 
Hamiltonian): 

Pp{W) = ip(W) exp[-pW\. (33) 

Then trivially extending the analysis of Ref. [4] we find 
that the inverse Fourier transform of the postselected 



work distribution function gives the return amplitude at 
a complex time: 

J dWPp [W) exp[iWt] = G(t + 0) (34) 

Thus amazingly the complex time return amplitude and 
hence the Fisher zeros analyzed in Fig. 1 of the main 
paper are measurable quantities at least in principle (note 
that the postselected temperature can be both positive 
and negative). 
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